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ABSTRACT 

As a step towards developing a new design 
philosophy, one that moves away from the 
traditional empirical approach used today in 
design towards a science-based design 
technology approach, a test series of 7 iso- 
tropic shells carried out by Arbocz and 
Babcock [1] at Caltech is used. It is shown how 
the hierarchical approach to buckling load 
calculations proposed by Arbocz et al [2] can 
be used to perform an approach often called 
“high fidelity analysis", where the uncertainties 
involved in a design are simulated by refined 
and accurate numerical methods. The Delft 
Interactive Shell DEsign COde (short, 
DISDECO) is employed for this hierarchical 
analysis to provide an accurate prediction of 
the critical buckling load of the given shell 
structure. This value is used later as a 
reference to establish the accuracy of the 
Level-3 buckling load predictions. As a final 
step in the hierarchical analysis approach, the 
critical buckling load and the estimated 
imperfection sensitivity of the shell are verified 
by conducting an analysis using a sufficiently 
refined finite element model with one of the 
current generation two-dimensional shell 
analysis codes with the advanced capabilities 
needed to represent both geometric and 
material nonlinearities. 
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INTRODUCTION 

It is generally agreed that, in order to make 
the development of the Advanced Space 
Transportation System a success and to 
achieve the very ambitious performance goals 
(like every generation of vehicles lOx safer and 
lOx cheaper than the previous one), one must 
make full and efficient use of the technical 
expertise accumulated in the past 50 years or 
so, and combine it with the tremendous 
computational power now available. It is 
obvious that with the strict weight constraints 
used in space applications these performance 
goals can only be achieved with an approach 
often called “high fidelity analysis”, where the 
uncertainties involved in a design are 
simulated by refined and accurate numerical 
models. In the end the use of “high fidelity” 
numerical simulation will also lead to overall 
cost reduction, since the analysis and design 
phase will be completed faster and only the 
reliability of the final configuration needs to be 
verified by structural testing. 

The light-weight shell structures used in 
aerospace applications are often buckling 
critical. The buckling load calculations are 
usually carried out by one of the many 
currently available finite element based 
computer codes [e.g., 3,4]. In order to reduce 
computer execution time, buckling analyses 
are often done using only the small 
displacement stiffness matrix K 0 . This 
approach is used, despite the fact that the 
“initial stability problem” so formulated can only 
give physically meaningful answers if the 
elastic solutions based on K 0 (at least 

approximately) are identically equal to zero [5]. 

When the qualitative nature of the expected 
behavior is completely unknown, the stability of 
the structure must be investigated using the full 
tangent stiffness matrix Kj in order to 

guarantee accurate and reliable buckling load 
and buckling mode predictions. In order to 
discover the load level at which K j ceases to 
be positive definite (that is, the load level when 
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buckling occurs), a step-by-step analysis 
procedure is needed. 

In addition, it is imperative (though often 
completely neglected), that at the beginning of 
any stability investigation, the accuracy of the 
discrete model used should be checked 
against available analytical or semi-analytical 
results. This step is part of a mandatory study 
needed in order to establish the dependence of 
the buckling load predictions on the mesh 
distribution used. Furthermore, as has been 
pointed out in the past by Byskov [6], if one 
carries out imperfection sensitivity 
investigations, which involve an extension of 
the solution into the postbuckling response 
region, further mesh refinement may be 
needed since the wavelength of the dominant 
large deformation pattern may often decrease 
significantly. 

Finally, whenever one is engaged in shell 
stability analysis it is especially important that 
one is aware of the possible detrimental effects 
of a whole series of factors, that have been 
investigated extensively in the late 1960s and 
the early 1970s. Thus for an accurate and 
reliable prediction of the critical buckling load 
of a real structure, one must account not only 
for the influence of initial imperfections [e.g., 
7,8] and of the boundary conditions [e.g. 9], but 
one must also consider the effects of stiffener 
and load eccentricity [e.g., 10] and the 

prebuckling deformations caused by the edge 
restraints [e.g. 11,12]. 

A test series of 7 isotropic shells carried 
out by Arbocz and Babcock [1] at Caltech is 
used to illustrate how such a hierarchical 
approach to buckling load calculations can be 
carried out. The platform for the multi-level 
computations, needed for an accurate 
prediction of the critical buckling loads and a 
reliable estimation of their imperfection 
sensitivity, is provided by DISDECO [13]. With 
this open ended, hierarchical, interactive 
computer code the user can access from his 
workstation a succession of programs of 
increasing complexity. 


SOLUTION OF THE BUCKLING PROBLEM 

In the following it will be shown that with the 
help of DISDECO, the Delft Interactive Shell 
DEsign COde, the shell designer can study the 
buckling behavior of a specified shell, calculate 
its critical buckling load quite accurately and 
make a reliable prediction of the expected 
degree of imperfection sensitivity of the critical 
buckling load. The proposed procedure 
consists of a hierarchical approach, where the 
analyst proceeds step-by-step from the simpler 
(Level-1) methods used by the early 


investigators to the more sophisticated 
analytical and numerical (Level-2 and Level-3) 
methods used presently. 

Level-1 Perfect Shell Buckling Analysis 

The geometric and material properties of 
the isotropic shell A-8 of Ref. [1] are listed in 
Table 1. 

Table 1. 

Geometric properties of Caltech isotropic shell A-8 


t total (= 

h) =0.00464 

in (=0.117856 

mm ) 

L 

= 8.0 

in (=203.2 

mm ) 

R 

= 4.0 

in (=101.6 

mm ) 

E 

= 15.2x10 6 

psi (= 1.0480x10 5 

N/mm 2 ) 

V 

= 0.3 




Assuming a perfect shell (W = 0) and the 
following membrane prebuckling state 


W (o) = hW v =hAi2- 
c 

F (o) = .E5il l 2 
cR 2 

where 


( 1 ) 


X = 


a N x Eh 

- 7. • °C£ - 77 ■ Ncf - a cf h 

a cf ™Cf c “ 


and c 


= >/3(1-v 2 ) 


then the nonlinear equations governing the 
prebuckting state are identically satisfied and 
the linearized stability equations reduce to a 
set of equations with constant coefficients. It 
has been shown in Ref. [14] that by assuming 
an asymmetric bifurcation mode of the form 


W (1 > = hsinmir—cosn— (2) 

L R 

where 

m = k = number of axial half waves 
n = i = number of circumferential full waves 

one can reduce the solution of the linearized 
stability equations to an algebraic eigenvalue 
problem. Notice that the eigenvalue X^^ 
depends on the wave numbers m and n. The 
critical load parameter X c is the lowest of all 
possible eigenvalues. Thus finding X c 
involves a search over the integer valued wave 
numbers m and n. Using the Level-1 
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computational module AXBIF [14] a search 
over integer valued axial half-wave numbers m 
yielded the lowest eigenvalues listed in Table 2 
for the specified circumferential wave numbers 
n. 

Table 2. 

Buckling loads of the Caltech isotropic shell A-8 
Buckling load map for the perfect shell 
using AXBIF [14] (N c ^ = - 49.51517 Ib/in) 


n = 4 

X™ = 1 000 

(m = 34) 

n = 5 

X? = 1 .000 

(m -- 34) 

n = 6 

Xc = 1 .000 

(m = 34) 

n = 7 

x™ = i .ooo 

(m = 33) 

n = 8 

x™ = i .ooo 

(m = 33) 

n = 9 

Xc = 1 .000 

(m * 1) 


X™ = i .ooo 

(m « 33) 

n = 10 

X™ = 1 .000 

(rn = 33) 

n = 11 

Xc = i .ooo 

(m = 33) 

n = 12 

x™= i.ooo 

(m = 32) 

n = 13 

Xc = 1008 

(rn = 2) 


x™ = i .ooo 

(m = 32) 

n = 14 

Xc = 1.000 

(m = 32) 

n = 15 

X^= 1.001 

(m = 3) 


x™ = i .ooo 

(m = 31) 

As is 

known the 

lowest eigenvalue 


repeated eigenvalue with a high degree of 
multiplicity. The family of modes belonging to 
the lowest eigenvalue X, c = 1.000 lie on the so- 
called Koiter circle [7], 

To facilitate the interpretation of the 
numerical results obtained, DISDECO provides 
the user with various graphical interfaces. 

Thus the results of the search for the critical 
(lowest) buckling load X c can be displayed in 
a contour map as shown in Fig. 1. In order to 
provide a quick overview of the distribution of 
eigenvalues, the values displayed in the 
contour plot are normalized Thus in Fig. 1 the 
following normalized eigenvalues are plotted 

K|H1 

X[P=— ^ 

C -49.51517 

Notice that the critical buckling load can be 
calculated using a simple multiplication 


N? = ^N cf = 

1.000(—49.51517) = —49.51517 Ib/in 

Level-2 Perfect Shell Buckling Analysis 

To investigate the effects of edge constraint 
and of different boundary conditions on the 
critical buckling load of the perfect shell 
(W = 0) one has to switch to the Level-2 

computational module ANILISA [15]. In this 
module the axisymmetric prebuckling state is 
represented by 

w(°) = hW v + hw Q (x) 

(3) 

F<°) = ^[-Axy 2 + R 2 f 0 (x)] 
cR 2 

It has been shown in Ref. [15] that with these 
assumptions the prebuckling problem is 
reduced to the solution of a single fourth order 
ordinary differential equation with constant 
coefficients, which always admits exponential 
solutions. Closed form solutions for simply 
supported and clamped boundary conditions 
have been published in the literature [16]. 

For isotropic shells the linearized stability 
equations admit separable solutions of the 
form 

= hw,(x)cosn0 

(4) 

p0) = fi(x)cosn6 
c 

y 

where 0 = — . 

Using a generalization of Stodola's method 
[17] first published by Cohen [18] the resulting 
nonlinear eigenvalue problem is reduced to a 
sequence of linearized eigenvalue problems. 
The resulting ordinary differential equations are 
solved numerically by a technique known as 
“parallel shooting over N-intervals” [19]. Notice 
that by this approach the effect of edge 
restraint and the specific boundary conditions 
are satisfied rigorously. To find the critical load 
parameter X c an n-search must be carried 
out, whereby one must be careful to find not a 
local minimum but the absolute minimum. As 
can be seen from the results presented in 
Table 3 the n-search using membrane 
prebuckling and a rigorous satisfaction of SS-3 
(N x =v = w = M x =0) boundary conditions for 
the stability problem yields indeed the repeated 
eigenvalues —1.0 located on the “Koiter" 
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circle. The slight variations in the eigenvalues 
are due to the fact that both the axial half-wave 
numbers m and the circumferential full wave 
numbers n are integers. 

Table 3 

Buckling loads of the Caltech isotropic shell A-8 
(N^ =—49.51517 Ib/in) 

Buckling load map for the perfect shell using 
ANILISA [1 5] (B.C. N x = v = w = M x = 0) 


Prebuckling: 

Membrane 

Nonlinear 

n = 8 

X^= 1.10977(m=1) 

Xc'= 0.951 155 


X™ = 1 .00001 (m=33) 

X^ =0.966181* 

n = 9 

X™ = 1.00004(m=1) 

Xc =0.934990 


CO 

CO 

II 

E 

o 

o 

o 

o 

o 

T" 

1) 

Xc = 0.957035* 

o 

ii 

c 

X™ = 1.08160(m=1) 

Xc =0.940209 


X™= 1.0001 1(m=33) 

X2'= 0.946124* 

n = 11 

X{P= 1.12539(m=2) 

Xp =0.931492* 


X[P= 1.00035(m=32) 

XC‘= 0.939878 

n = 12 

X[P= 1.01 481 (m=2) 

X2’= 0.915200* 


Xp = 1 .00004(m=32) 

Xc =0.931464 

n = 13 

X™ = 1 ,00836(m=2) 

Xc= 0.908528* 


Xc = 1 .00004(m=32) 

Xc =0.918482 

n = 14 

X{T= 1.04064{m=3) 

Xc= 0.903256 


X™ = 1.00030(m=31) 

X2‘= 0.906963* 

n = 15 

Xc = 1.00060(m=3) 

Xc =0.891368 


X™ = 1.00000(m=31) 

Xc =0.898837* 

n = 16 

X?= 1.02041 (m=3) 

X2'= 0.885185 


X™= 1.00023(m=31) 

Xc =0.886606* 

n = 17 

X[T = 1.00086(m=4) 

Xc= 0.875881* 


X™ = 1.00001 (m=30) 

Xc= 0.878312 

n = 18 

X™ = 1.01300(m=4) 

X2'= 0.868406* 


XjTU 1.00027(m=30) 

Xc= 0.869200 

n = 19 

X™= 1.0001 4(m=5) 

Xc =0.861087 


X™= 1.00001 (m=29) 

X^= 0.861626* 

o 

C\J 

II 

C 

X™ = 1.00164(m=6) 

Xg 1 = 0.854891 


X^= 1.00006(m=28) 

X2'= 0.855033* 


Table 3 - continuation 


Prebuckling: 

Membrane 

CM 

II 

c 

X[P= 1.00378(m=7) 
X™= 1.00023(m=27) 

n « 22 

X[P= 1.001 83(m=7) 
X™ = 1.0001 4(m=27) 

n = 23 

X™ = 1.001 14(m=8) 
X{P= 1.0001 2(m=26) 

CM 

II 

C 

Xc= 1.00071 (m=10) 
X™ = 1.0001 1(m=24) 

n = 25 

X™ = 1.00000(m=11) 
X™= 1.00000(m=23) 

n = 26 

X|P= 1.00002(m=13) 
XjP= 1 . 00001 (m=21) 

n = 27 

X™ = 1 ,00027(m=1 7) 
Xg 1 = 1.00032(m=18) 


Nonlinear 

Xq = 0.849839* 
Xc l = 0.849927 
Xg' = 0.846336* 
X2'= 0.846344 
X2'= 0.844502 
= 0.844511* 
X^ 1 =0.844480* 
0.844481 
X2'= 0.846279* 
Xc = 0.846280 
X^ 1 =0.849833* 
Xc =0.849833 
X2'= 0.854990* 
X2'= 0.854990 


* anti-sym at x = L/2. 

The most accurate Level-2 solutions are 
obtained when one employs a rigorous 
nonlinear prebuckling analysis. As can be seen 
from the results listed in Table 3, for this 
particular shell the critical buckling loads with 
nonlinear prebuckling are always lower than 
the corresponding results obtained using a 
membrane prebuckling analysis. Notice that 
the critical load N c can be calculated easily by 
multiplying the lowest eigenvalue Xq by the 
normalizing factor Nc/ = -49.51517 Ib/in 
yielding 


Nc = ^Nctf = -41 .81 5 Ib/in (n = 24) 

In fig. 2 the critical buckling modes using 
membrane and rigorous nonlinear prebuckling 
are depicted. Notice that the solutions with 
nonlinear prebuckling differ significantly from 
the ones obtained using membrane 
prebuckling, especially at n = 24 where one 
observes a typical edge buckling type 
behavior. 

Level-3 Perfect Shell Buckling Analysis 
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To verify the earlier predictions the finite 
difference version [20] of the well known shell 




analysis code STAGS [21] will be used. In 
order to be able to represent the measured 
initial imperfections accurately, the whole shell 
will be modeled. 

Initially a convergence study must be 
carried out in order to establish the mesh size 
needed for accurate modeling of the buckling 
behavior of the shell in question. For this 
purpose the asymmetric bifurcation from a 
nonlinear prebuckling path option was used, 
whereby the earlier results obtained with the 
Level-2 module ANILISA listed in Table 3 
serve as a reference. 

In the convergence study, at first, for a fixed 
number of mesh points in the axial direction 
(NR = 161) the number of mesh points in the 
circumferential direction (NC) was increased 
until the bifurcation load approached a 
horizontal tangent. As can be seen from Fig. 3 
the results converge to a limiting value from 
below at about NC = 261. Next, for a fixed 
number of mesh points in the circumferential 
direction (NC = 161) the number of rows (NR) 
was varied. This time convergence is from 
above and as can be seen from Fig. 3 the 
horizontal tangent is reached at about NR = 
261. 

To illustrate the difference between using 
coarser meshes to speed up the computations 
and the finer meshes which produce (nearly) 
converged solutions, initially the first three 
eigenvalues and eigenvectors are computed. 

Using a mesh of 161 rows and 181 columns 
(a model with 89488 D.O.F.’s and a maximum 
semi-bandwidth of 713) and SS-3 boundary 
conditions (N x = -N 0 ,v = W = M x =0) the 
following 3 lowest eigenvalues were obtained: 

4 1 ) = 0.833361 (n = 26) -4 

= 4 1 ) n c ^ = “ 41 - 264 lb/in 


= 0.844993 (n = 24) -» 

= " 41 840 lb/in 

4 2 > = 0.845006 (n = 24) 

=- 41 -841 lb/in 

~ 0.8451 00 (n = 25) -> 

Nc 3) = 4 3 H* = - 41 845 lb/in 

Details of the critical mode are displayed in 
Figures 4-5 for the coarser 161x181 model and 
in Figures 6-7 for the finer 261-261 model. 
Notice that using the more refined 261x261 
model the sequence of the 3 lowest buckling 
loads and the corresponding buckling modes 
agree closely with the predictions obtained with 
the Level-2 module ANILISA (see also Table 3) 


IMPERFECTION SENSITIVITY STUDY 

That initial imperfections may decrease the 
load carrying capacity of thin-walled shell 
structures is by now widely known and 

accepted. However, in order to calculate the 
effect of initial imperfections one must know 
their shape and amplitude, an information that 
is rarely available. 

In the absence of initial imperfection 

measurements, as a first step one must 
establish whether a given shell-loading 

combination is imperfection sensitive, and if 
the answer is positive to estimate how 
damaging certain characteristic imperfection 

shapes are. 


A*. 2 ) = 0.833374 (n = 26) -» 

n c 2 > = ^c 2 ^ N c^ =-41.265 lb/in 

= 0.834238 (n = 25) -* 

N< 3 > = = -41 .307 lb/in 


Single Axisymmetric Imperfection 

Based on Koiter's pioneering work on the 
effect of initial imperfections [7,22] the simplest 
imperfection model consists of a single 
axisymmetric imperfection 


W = h^j cosin' 


(5) 


Next these computations were repeated using 
a mesh of 261 rows and 261 columns (a model 
with 207508 D.O.F.’s and a maximum semi- 
bandwidth of 1037) and the same SS-3 
boundary conditions yielding the following 
results 


where i is an integer denoting the number of 
half-waves in the axial direction and si is the 
amplitude of the axisymmetric imperfection 
normalized by the shell wall-thickness h. 

If one assumes that both the axial load and the 
boundary conditions are independent of the 
circumferential coordinate, then the prebuck- 
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ling solution will also be axisymmetric, a fact 
that simplifies the solution considerably. 


Sj = 2m if i = 2m ! 

= 0 otherwise 


Level-1 Analysis of Axisymmetric 

Imperfection 

Neglecting the effect of the prebuckling 
boundary conditions the nonlinear equations 
governing the prebuckling state admit the 
following axisymmetric solutions 


W<°) - hW v + hw 0 (x) 


where 


( 6 ) 


Wo(x) = — - ^ h^ 1 cosing 


fo(x) = 


X c - -X 


X (1 + B2iap) E h 3 

k c\ 2apA22 c 


cos Its- 


1 2 _* (1 + B2iaf) 2 

X o\ =2{a, 2 Dn+ — 4 } 

2 af A 22 


(7) 


Notice that the linearized stability equations 
become now a set of equations with variable 
coefficients. The reduced wave number otj 

* 

and the normalized stiffness coefficients A 22 » 
♦ * 

B 21 and Di i are all listed in Ref. [14]. 

It has been shown in Ref. [14] that by 
assuming an asymmetric bifurcation mode of 
the form 

= hsinrrm— cos n (8) 

L R 

where 

m = k = number of axial half waves 
n = e = number of circumferential full waves 


5j m =1 if i = m 
= 0 otherwise 

and the constants C-|,C2,... are listed in 
Appendix C of Ref. [14]. 

Here it must be remembered that one will 
only get any noticeable degrading influence of 

the assumed axisymmetric imperfection if ^ 

is negative and if the coupling condition i=2m 
is satisfied. The physical explanation for this 
can be found in Koiter’s 1963 paper [22]. 
Furthermore, in order to obtain the smallest 
real root of Eq. (9), for a given axisymmetric 
imperfection ^ an n-search must be carried 

out. It should also be noticed that the terms 
involving the Kronecker delta Sj = 2m are 

linear in ^ , and thus they dominate the buck- 
ling behavior of the shell with axisymmetric 
imperfection. 

Assuming that the most likely axisymmetric 
imperfection of the silver pointed wax mandrel 
used to electroplate the isotropic shell A-8 is 
given by 


W = h^-j cos 271-^- (10) 

the Level-1 OISDECO computational module 
AXBIF generated the solid curve shown in Fig. 
8. Notice that the curve is normalized by 
N c ^ = -49.51517 Ib/in, the classical buckling 
load of an isotropic shell, obtained around 
1910 by Lorenz [23], Timoshenko [24] and 
Southwell [25]. Notice also that an initial 
imperfection amplitude equal to the wall 
thickness of the shell (5i = -1 .0) generates a 
“knockdown factor” of Xc = 0.706 , resulting in 
the following buckling load 


a Galerkin type approximate solution yields for 
the eigenvalues (read, buckling loads) X of the 
problem a characteristic equation in the form of 
a cubic polynomial 

^ - (^-mnt + 2^-cj - Ci4-j8j=2m)^ (9) 

+ {2A. mrn + A, C j + (62 - C-|)^i8j = 2m}^Cjk 

- {^mnt + ^2^i=2m + (C3 - 648^)^ }X 2 = 0 


N c — XcNc* — 

0.706(—49.51517) = -34.958 Ib/in 

However, if one assumes an imperfection in 
the form of the classic axisymmetric buckling 
mode of Koiter [7] 

W = cosi^n^ 
where 


where 
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( 12 ) 


L [2c 
' c ' n VRh 


c = ^3(1- v 2 ) 


then for the isotropic shell A-8 i c ^ - 34.0 , and 
the Level-1 DISDECO computational module 
AXBIF found that the minimum buckling load 

occurred when m = ^cf 1° r different values of 

n (which depended on the value of ^ ). These 

results are shown in Fig. 9. 

For this case an initial imperfection amplitude 
equal to the wall thickness of the shell 
= -1 .0) generates a “knockdown factor” of 
Xc = 0.196 resulting in the following rather low 
buckling load 


=hwi(x)cosn6 

p( 1 )=§5^f 1 (x)cosne 
c 

y 

where 0 = r-. 

R 

Solution proceeds as outlined on Ref. [26]. 
Using an updated version of the Level-2 
computational module MAN I LISA [27] and SS- 
3 (N x = v = w = M x =0) boundary conditions 
one obtains the results presented in Table 4. 
Notice that a rigorous nonlinear prebuckling 
analysis was used and an n-search was 
carried out for each specified axisymmetric 

imperfection amplitude \ -j . 


Nc = X^ = 

0.196(-49.51517) = -9.705 Ib/in 

The last two buckling loads clearly 
demonstrate the importance of the shape of 
the initial imperfections involved and the role of 
the corresponding eigenvalues. Using the long 
wave axisymmetric imperfection (i = 2) the 
corresponding normalized axisymmetric 
buckling load (eigenvalue)of the perfect shell is 
Xc. =144.321 (a very high value) yielding a 

low imperfection sensitivity. For the short wave 
imperfection (i = 34) the corresponding 
normalized axisymmetric buckling load 
(eigenvalue) of the perfect shell is 
Xc. =1.000. Thus it lies on Koiter’s circle 

yielding a high imperfection sensitivity. 

Level-2 Analysis of Axisymmetric 

Imperfection 

Since the external loading, the boundary 
conditions and the assumed initial imperfection 
are axisymmetric, therefore the prebuckling 
solution will also be axisymmetric. It has been 
shown in Ref. [26] that by assuming 

w(°) = hWv + hw D (x) 

( 11 ) 

F < 0) =^-{-ix y 2 +r 2 , oW> 

the solution of the nonlinear partial differential 
equations governing the prebuckling state can 
be reduced to the solution of a single fourth 
order ordinary differential equation with 
constant coefficients, which can be solved 
routinely. 

For isotropic shells the resulting linearized 
stability equations admit separable solutions of 
the form 


The values of Table 4 are plotted as the 
dashed curve in Fig. 9. A comparison of the 
results obtained via the Level-1 module AXBIF 
(solid curve) and the Level-2 module 
MANILISA (dashed curve) shows that also in 
the case of axisymmetric imperfections a 
rigorous pre-buckling analysis should be used. 
Especially for very small initial imperfection 

amplitudes (|^-|| < 0.1) the Level-1 predictions 

are inaccurate and overestimate the critical 
buckling load. Notice further that both curves 
have been normalized by classical buckling 
load N c ^ = -49.51517 Ib/in. This way the effect 
of using a rigorous prebuckling analysis 
becomes easily discernible. 

Table 4 

Buckling loads of the Caltech isotropic shell A-8 
(N c ^ = —49.51517 Ib/in) 

Short wave axisymmetric imperfection 


W/h = ^cos34tc^ using MANILISA [27] 


(B.C. : N 

x =v = w = M x = 

0) 


Si 

X? 

Si 

^nl 

0 . 

0.844480 (n=24) 

-0.1 

0.549168 (n=17) 

-0.01 

0.847663 (n=21 ) 

-0.2 

0.420091 (n=17) 

-0.02 

0.786079 (n=21 ) 

-0.3 

0.354645 (n=19) 

-0.03 

0.740449 (n=21) 

-0.4 

0.305484 (n=18) 

-0.04 

0.703421 (n=21 ) 

- 0.5 

0.271507 (n=17) 

-0.05 

0.671958 (n=21) 

-0.6 

0.247002 (n=16) 

-0.06 

0.644481 (n=21) 

-0.7 

0.229213 (n=1 5) 

- 0.07 

0.620047 (n=21) 

-0.8 

0.214064 (n=20) 

- 0.08 

0.598039 (n=21) 

-0.9 

0.200103 (n=19) 

-0.09 

0.578030 (n=21) 

- 1.0 

0.189011 (n=18) 
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Level-3 Analysis of Axisymmetrlc 

Imperfection 

Recalling that since both the axial load and 
the boundary conditions are independent of the 
circumferential coordinate, therefore the 
prebuckling solution will also be axisymmetric, 
one can use once again the asymmetric 
bifurcation from a nonlinear prebuckling path 
option. By modeling the full shell the code can 
choose itself the critical number of full waves in 
the circumferential direction. No n-search must 
be carried out. Using a uniformly spaced mesh 
of 161 rows and 181 columns and the user 
written subroutine option WIMP to introduce 
the following axisymmetric imperfection 

W = 10 h cos 2k — 

L _ 

(remember STAGS defines W positive 
outward) the following critical bifurcation load 
was found 

N c = -32.582 Ib/in 

As can be seen from Fig. 10 the critical 
buckling mode has n = 9 full waves in the 
circumferential direction and a single half wave 
in the axial direction. The nondimensional 
bifurcation load of the shell with axisymmetric 
imperfection is for ^ = 1.0 


circumferential full waves, respectively, 
instability occurs at the limit point of the 
prebuckling state in the generalized load- 
deformation space. Assuming that the 
eigenvalue problem for the critical (lowest) 
buckling load A c will yield a unique asym- 
metric buckling mode W^ f then for an 
imperfect shell (£ 2 * 0) the shape of the 
generalized load-deflection curve in the vicinity 
of the bifurcation point A = A c is given by the 
following asymptotic expansion 

(A - A c )£ = A c a£ 2 + A c b£ 3 + ... 

(14) 

-A c a$ 2 — (A — A c )05 2+ O(i£ 2 ) 

Expressions for the postbuckling coefficients 
“a" and “b" and the imperfection forms factors 
“a” and “|T are derived in References [28,29]. 
If the limit point is close to the bifurcation point, 
then the maximum load A s that the structure 
can carry prior to buckling can be evaluated 
from Eq. (14) by maximizing A with respect to 
q . For cases where the first postbuckling 
coefficient “a” is zero, this analysis yields the 
modified Koiter formula [29] 

(1-Ps) 3/2 = |V-3a 2 b[1-£(1-p s )]|f 2 |(15) 


X _ Nc _ -32.582 
C N c ^ -49.51517 


0.65802 


Notice also that the Level-2 MANILISA 
prediction ()$ = 0.683, n = 9) agrees closely 
with the Level-3 STAGS-A prediction 
(A,c* = 0.658, n = 9) . The slight difference is 
partly due to the fact that MANILISA uses the 
Donnell type nonlinear shell equations, 
whereas STAGS-A employs the higher order 
Marlowe-FICigge equations. 

Single Asymmetric Imperfection 


where p s = A s /A c . 

Notice that, if the second postbuckling 
coefficient “b M is positive, Eq. (15) has no real 
solutions. Thus the buckling load of the 
specified shell-loading combination is not 
sensitive to small asymmetric initial 
imperfections of the shape given by Eq. (13). 
If, however, the second postbuckling 
coefficient “b" is negative, the equilibrium load 
A decreases following buckling and the 
buckling load of the real structure A s is 
sensitive to the asymmetric initial imperfection 
specified by Eq. (13). 


The effect of a single asymmetric initial 
imperfection can be investigated either by 
solving the full nonlinear response problem or 
by employing the well known Lyapunov- 
Schmidt-Koiter [7] reduction technique. When 
investigating the degrading effect of a single 
mode asymmetric imperfection 

W = t^p sinrrm— cosn-^- (13) 

L R 

where m and n are integers denoting the 
number of axial half-waves and the number of 


Level-1 Analysis of Asymmetric 

Imperfection 

For the isotropic shell under investigation, 
as can be seen from the partial results listed in 
Table 2, there are many identical eigenvalues. 
Hence, strictly speaking, the proposed form of 
the perturbation expansion given by Eqs. (14) 
is not applicable, since the nonlinear 
interaction between the many nearly 
simultaneous eigenmodes is not accounted for. 
Thus the following results, where one 
considers the eigenfunctions corresponding to 
certain critical eigenvalues chosen one at the 
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time, can at best give an indication as to the 
severity of the expected imperfection 

sensitivity. 

Initially, we assume a long wave 

asymmetric imperfection affine to one of the 
critical buckling modes of the perfect Caltech 
isotropic shell A-8, computed by the Level-1 
computational module AXBIF (see also Table 
3) 

w =h^2 s ' n7t r cos9 o‘ ( 16 ) 

L R 

Next using the Level-1 computational module 
BFACT to carry out the initial postbuckling 
analysis yields the following results 

X c = X™ = 1 .00004 (m = 1, n = 9) 

b = -0.010636 a = (1 = 1.0 

Substituting these values into Eq. (15), one 
can plot the degrading effect of an asymmetric 
imperfection of the shape given by Eq. (16) as 
a function of its amplitude £2 • As can b® seen 
from Fig. 11 an initial imperfection amplitude 
equal to the wall thickness of the shell 
(£ 2 = t0) generates a “knockdown factor” of 
X s = 0.679 , resulting in the following low 
buckling load 

N s = ^sNc^ 

0.679 (-49.51 517) = -33.621 Ib/in 

Notice that the imperfection form factors 
"a" and "P“ are identical equal to 1.0 because 
BFACT uses membrane prebuckling to 
calculate the necessary first and second order 
fields and the assumed asymmetric 
imperfection shape of Eq. (16) is affine to the 
buckling mode. If, however, one assumes the 
following short-wave asymmetric imperfection 

W - hlosin17jt^-cos27^- 
L R 

which lies also on the Koiter circle, then the 
initial postbuckling analysis, done once again 
by the Level-1 computational module BFACT, 
yields the following results 

Xc=Xc = 1 .00025 (m = 1 7, n = 27) 

b = -220.84 a = P = 1 .0 

Using these values in Eq. (1 5) one obtains Fig. 
12. As can be seen, an initial imperfection 


amplitude equal to the shell wall thickness 
(^2 = 1 .0) this time generates a "knockdown 
factor” of X s = 0.0249 , resulting in an 
unrealistically low buckling load of 

N s = X«Ncf = 

0 . 0249 {— 49.51517) = -1 .233 Ib/in 

Here one obviously cannot rely on the 
prediction of the Level- 1 computational module 
BFACT, which uses membrane prebuckling 
analysis. In such a situation, as has been 
shown by Hutchinson and Frauenthal [30], one 
must repeat the postbuckling analysis using a 
rigorous prebuckling solution which includes 
the nonlinear effects due to the edge restraints. 

Level-2 Analysis of Asymmetric 

Imperfection 

To investigate the effects of edge-constraint 
and/or different boundary conditions on the 
imperfection sensitivity of the critical buckling 
load of the Caltech isotropic shell A-8 one has 
to switch to the Level-2 module ANILISA [15] 
and run its postbuckling analysis option. In this 
module, as described earlier, the axisymmetric 
prebuckling state is represented by Eqs. (3), 
the buckling modes by Eqs. (4) and the 
postbuckling state by 

W< 2 > = h[w a (x) + wp(x)cosn0] 

(18) 

2 

p(2) =I^L[f a (x) + fp(x)cos2n0] 

where 0 = y/R. Details of the computational 

procedures used are reported in Refs. [15,26]. 

Initially, let us assume that the specified 
asymmetric imperfection is given by Eq. (16) 

W = h^ 2 s i njt ^ cos 98 

where 0 = y/R . Running ANILISA with rigo- 
rous prebuckling and SS-3 boundary 
conditions (N x = -No, v = w = M x = 0) yields 
the following results 

Xc = Ag 1 = 0.934990 (n = 9) 

b = -0.016202; a =0.20830; P = -0.38232 

Using Eq. (15) to plot the degrading effect 
of the asymmetric imperfection specified by 

Eq. (16) as a function of its amplitude %2 one 
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obtains the results displayed in Fig. 13 as a 
solid line. 

Here one must remember that Koiter’s 
Sensitivity Theory is asymptotically exact, that 
is, it yields accurate predictions for sufficiently 
small imperfections, whereby what is 
sufficiently small may vary from case to case. 
Also, Eq. (15) was obtained by using the 
perturbation expansion given by Eq. (14), 

where terms of order (££) are neglected. As 
can be seen from the dotted curve plotted in 
Fig. 13, by using more advanced 
computational modules such as COLLAPSE 
[31], where a full nonlinear solution is used and 

-p 

terms up to and including order (££ ) are 
kept, one obtains slightly lower predictions. 

Notice that up to about £2 = 01 the 
asymptotic predictions from ANILISA and the 
nonlinear results of COLLAPSE agree very 
closely. Thus one can say that in this case the 
range of validity of the asymptotic solution is 

o>i 2 so.i. 

For an initial imperfection amplitude equal 
to the wall thickness of the shell (^ 2 = 1 .0) the 
asymptotic solution yields a “knockdown factor” 
of X s = 0.806 , which corresponds to a 
buckling load of 

Ns=^sNc* = 

0.806(-49.51517) = -39.909 Ib/in 

For the same initial imperfection amplitude 
of ?2 = 1 0 nonlinear solution using 

COLLAPSE yields a “knockdown factor'’ of 
X s = 0.68 , which corresponds to a buckling 
load of 

Ns — X s Nc^ = 

0.68{— 49.51517) = -33.670 Ib/in 

Let us next consider the short wave 
asymmetric imperfection given by 

W = h ^2 sin177i^cos24G (19) 

where 0 = y / R . The rigorous asymptotic 
solution obtained by running ANILISA with SS- 
3 boundary conditions (N x =-N 0 , 
v = w - M x =0) yields now the following 
results 

x c = Xg 1 = 0.844480 (n = 24) 


b = -0.50329; a = 0.201 13; 0 = -0.012022 

Using these values with Eq. (15) yields the 
solid curve plotted in Fig. 14. Rerunning the 
same case with the nonlinear module 
COLLAPSE yields the dotted curve plotted in 
Fig. 14. 

Notice that in this case the asymptotic 
predictions from ANILISA and the nonlinear 
results from COLLAPSE agree closely up to 
about ^ 2 =0.03. Thus the range of validity of 

the asymptotic solution is only 0 £ ^ < 0.03 . 

For an initial imperfection amplitude of 
^2 = 1 .0 the asymptotic solution yields a 
“knockdown factor" of k s = 0.400 , which 
corresponds to a buckling load of 

N s - * S N C ^ = 

0.400(— 49.51517) = -19.806 Ib/in 

For the same initial imperfection amplitude of 
^2 = 1 .0 the nonlinear solution yields a 
“knockdown factor" of X s = 0.30 , which 
corresponds to a buckling load of 

N s = X s Hq£ - 

0.30( -49.51 517) = -14.855 Ib/in 

For an initial imperfection amplitude equal 
to the wall thickness of the shell (^ 2 = 1 .0) the 
asymptotic solution yields a “knockdown factor" 
of k s =0.806, which corresponds to a 
buckling load of 

N s =*sNc* = 

0.806(-49.51517) = -39.909 Ib/in 

MEASURED INITIAL IMPERFECTIONS 

In order to apply the theory of imperfection 
sensitivity with confidence, one must know the 
type of imperfections that occur in practice. In 
1969 Arbocz and Babcock [1] published the 
results of buckling experiments where, for the 
first time, the actual initial imperfections and 
the prebuckling growth of the midsurface of 
electroplated isotropic shells were carefully 
measured and recorded by means of an 
automated scanning mechanism. 

Midsurface Initial Imperfections 

As can be seen from Fig. 15, the measured 
initial midsurface imperfections of shell A-8 [1] 
show a rather general distribution dominated 
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by an n=2 mode. One can use the following 
double Fourier series 

^ = £w jo cosijc£ + £w 0< cos^ (20) 

" i=i r=i 

1 9 

+ £w^sin^ + XX^kfCO sk ^ cos ^ 

r=l M=1 

115 _ xv 

+ XZ WkfC 0 Sk 7 l L S ' n ^R 

M*1 


MANILISA [27] Xg 1 = 0.266 (n = 27) 
P" 1 = awRXg'Nc/ = -331 .024 lbs 

Turning now to the long wave asymmetric 

imperfection W = h ^2 sinn— cos9^- 
if \ z = Agsy = 0.6882 then from 

COLLAPSE [31] Xg 1 = 0.740 (n = 9) 
P" 1 = 2jiRXg*N c ^ = -920.894 lbs 


to represent the measured initial imperfections 
accurately, where Fourier coefficients with 
absolute values less than 0.001 are neglected. 

A convenient measure of the size of the 
initial imperfections is their root-mean-square 
(RMS) value. By definition 


For the short wave asymmetric 

x y 

imperfection W = h^g sin 1771— cos 24— 

if ^2 = A asy = 0.6882 then from 

COLLAPSE [31 ] X" 1 = 0.324 (n = 24) 


2tiRL 


a 2 _ 

^rms “ 


2nRL 


J J[W(x,y] 2 dxdy (21) 


o o 


Thus 




) 


M 
2 , a 2 


= A axi + Aas y 

For the A-8 shell the RMS values are 


A axi 


= 0.5077 


A aS y = 0.6882 


Using these measured RMS values with the 
previously discussed initial imperfection 
models yields the following buckling load 
predictions: 

For the long wave 

^ 

axisymmetric imperfection W =h^iCOs2n— 

if ^ = A ax j = -0.5077 then from 

MANILISA [27] Xg 1 = 0.824 (n = 9) 
P£' = 27iRXg , N c ; = — 1 025.428 lbs 


For the short wave 

x 

axisymmetric imperfection W = cos347t— 

if = A ax j = -0.5077 then from 


p£ ! = 2^X2^ = -403.202 lbs 


Comparing these buckling load predictions with 
the experimental buckling load of 

p ex p = -782.7 ,bs it is clear that the long wave 

imperfection models yield an upper bound, 
whereas the short wave imperfection models 
yield a lower bound. 

For a more accurate estimate of the 
buckling load one must use the measured 
initial imperfections in codes like STAGS 
[20,21] to carry out 2-dimensional nonlinear 
collapse analysis. Employing a user written 
subroutine WIMP to input the double Fourier 
series of Eq. (20), the 161x181 STAGS model 
yielded with SS-3 boundary condition 
(N x = -No, v = w = M x = 0) a collapse load of 

p s _ _900.103 ,bs , whereas the same model 
with C-4 boundary condition (u = u Q , 
v = w = w, x = 0) yielded a collapse load of 

p s = -976.467 ,bs . Both of these values are 
significantly higher than the experimental 

Ihe 

buckling load of P exp = -782.7 D . 

Boundary Initial Imperfections 

In References [32,33] it was shown that 
boundary imperfections can have a significant 
degrading effect on the buckling loads of 
axially compressed cylindrical shells. The 
flatness of the end ring, attached to the load- 
cell of the Caltech test set-up shown in Fig. 16, 
was measured about 20 years after it was 
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used to test shell A-8. As can be seen in Fig. 
17 there is a very large amplitude waviness of 
the end support, with maximum deviations of 
about + 3 wall thicknesses. 

The measured boundary imperfections are 
decomposed in a 1 -dimensional Fourier series 


u 0 =Ub(y)= (23) 

l bh | \ a 0 + £ (a n cos n ^ + b n sin n J 

This series then is used in STAGS [20] to 
model the effect of boundary initial 
imperfections using a modified C-4 boundary 
condition 

u = u b (y) , v = w = w, x = 0 (24) 

and taking advantage of the dual loading 
systems provided by the program. 

Because the grooves of the end rings were 
filled with liquid Cerrolow at the time shell A-8 
was installed in the test set-up, it is to be 
expected that after cooling off the hardened 
liquid metal has filled-in all the gaps. The effect 
of this unknown end support was modeled by 
varying the amplitude of the boundary 
imperfection between 0 and 0.1 . As can be 

seen from the results displayed in Fig. 18, 
including both the measured initial midsurface 
imperfections of Eq. (20) and the measured 
boundary imperfections of Eq. (23) with an 
amplitude of ^ - 0.035 , the calculated 

collapse load agrees well with the experimental 
buckling load. 


DISCUSSION OF THE RESULTS 

When comparing and analyzing the results 
obtained sofar it is important to keep in mind 
that all Level-1 and Level-2 solutions are 
based on approximate representations of the 
unknown functions. As pointed out in the 
previous sections Level-1 solutions use a 
single term double Fourier series 
approximation to reduce the solution of the 
stability problem, formulated in terms of partial 
differential equations, to algebraic eigenvalue 
problems. The effect of edge restraint is 
neglected (one uses a membrane prebuckling 
solution) and the assumed field functions 
satisfy approximately SS-3 (N x =-N 0 , 

v = w » M x = 0) boundary conditions. 

Level-2 solutions eliminate the y- 
dependence by a truncated Fourier 


decomposition in the circumferential direction. 
The resulting system of nonlinear ordinary 
differential equations are solved numerically, 
whereby both the specified boundary 
conditions and the effect of edge restraint are 
rigorously satisfied. Thus by this approach the 
only approximation is that one represents the 
variation of the solution in the circumferential 
direction by a single harmonic with n full 
waves, whereby an n-search is used to 
establish which wave number is the critical 
one. The Level-2 module ANILISA can also be 
used to investigate the effect of using different 
boundary conditions. In Table 5, the results for 
four different boundary conditions are 
presented. As expected the fully clamped C4 
boundary conditions has the highest critical 
buckling load. The increase in load carrying 
capacity with respect to the weaker SS3 
boundary conditions is about the same as for 
an isotropic shell of similar characteristic 
dimensions (same L/R and R/t ratios) of Ref. 
[34], 

The Level-3 solutions are based either on a 
2-dimensional finite difference or finite element 
formulation. In both cases, if one uses the 
appropriate meshes, one can obtain rigorous 
solutions where all nonlinear effects are 
properly accounted for. The only real problem 
with Level-3 type solutions is that for each 
problem one must establish the appropriate 
mesh size. Coarse meshes yield inaccurate 
solutions. What is coarse depends on the 
particular problem under investigation. Thus, 
for a general nonlinear solution a convergence 
study must always be carried out. 

Using a hierarchical simulation platform 
such as DISDECO (Delft Interactive Shell 
DEsign Code), where the analyst has at his 
disposal computational modules of different 
level of sophistication, such a convergence 
study can be carried out relatively quickly and 
accurately. In Table 5 a summary of the results 
obtained in this study is presented using 
normalized variables. 

Table 5 

Summary of buckling load calculations of the 
Caltech isotropic shell A-8 

(N c ^ = —49.51517 IbAin) 

Perfect shell analysis 

Level-1 Membrane Prebuckling - 

Repeated eigenvalues with high multiplicity 

X? = 1 .000 P c = -1 244.5 lbs 
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Level-2 Nonlinear Prebuckling 

SS-3 Xg 1 = 0.844480(n = 24) 

-» P c =-1050.9 lbs 
SS-4 Xg 1 = 0.866974(n = 25) 

P c =-1078.9 lbs 
C-3 Xg 1 = 0.910305(n = 24) 

P c =-1132.8 lbs 
C-4 Xg 1 = 0.926871(n = 25) 

-> P c =-1153.4 lbs 

Level-3 SS-3 B.C. 

161x181 Model 

Xg 1 = 0.833361(n = 26) 

-► P c =-1037.1 lbs 

261x261 Model 

Xg 1 = 0.844993(n = 24) 

-> P c =-1051.6 lbs 

Imperfect Shell Analysis 

Long wave axisymmetric imperfection - 

5l = Aaxi = -0-5077 

Level-2 Xg 1 = 0.824(n = 9) 

-> P c =-1025.4 lbs 

Short wave axisymmetric imperfection - 

?1=A axi = -0.5077 

Level-2 Xg 1 = 0.266(n = 27) 

-► P c = -331 ,0 lbs 

Long wave asymmetric imperfection - 

\z - A asy = 0-6882 

Level-2 Xg 1 = 0.740(n = 9) 

-> P s = -920.9 lbs 

Short wave asymmetric imperfection - 

?2 ~ A a sy — 0.6882 

Level-2 Xg' = 0.324(n = 24) 

P s = — 403.2 lbs 

Measured midsurface initial imperfection - 

see Eq. (20) 

Level-3 SS-3 B.C. -> P s = -900.1 ,bs 


Level-3 C-4 B.C. -> P s = -976.5 ,bs 

Measured midsurface + boundary (c^ = 0.035) 
initial imperfection - see Eqs. (20) and (23) 

Level-3 Modified C-4 B.C. 

(see Eq. (24)) -» P s = -780.0 lbs 

Experimental Buckling Load 

-> P eX p = -782.7 lbs 

Notice that by including both the midsurface 
and the boundary initial imperfections one 
obtained a very good agreement indeed. 

The study of the buckling behavior of the 
isotropic shell A-8 was undertaken with the 
objective of establishing a rational analysis 
procedure for incorporation of initial 
imperfections into the design of buckling critical 
structures. For a shell of revolution under 
axially symmetric loading there are the 
following four items necessary for the 
application of the proposed procedure: 

a. Buckling analysis capability for shells of 
revolution under axially symmetric loading. 
(A computer code such as SRA [18] or 
BOSOR [35] will satisfy this requirement.) 

b. Imperfection sensitivity analysis capability 
for shells of revolution under axially 
symmetric loading. (A computer code such 
as SRA [18] or STASOR [36] will satisfy this 
requirement.) 

c. Geometrically nonlinear analysis capability 
for a shell of revolution, which can handle 
the asymmetry introduced by the general 
asymmetric initial imperfections. (The 
STAGS [20,21] computer code will satisfy 
this requirement.) 

d. A design imperfection. 

At the present time this is the most elusive 
item. 

The imperfection model introduced by 
Imbert [37], following an idea by Donnell and 
Wan [38], is one approach, but has yet to be 
extended to larger cylindrical shell structures 
and other shells of revolution. 


CONCLUSIONS 

By relying on a series of theoretical results of 
various degree of sophistication published in 
the literature, the hierarchical approach used in 
this paper has resulted in a series of buckling 
load predictions of increasing accuracy. It was 
shown that in order to be able to arrive at a 
reliable prediction of the critical buckling load 
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and to make an estimate of its imperfection 
sensitivity which can be used with confidence, 
one must proceed step by step from simple to 
more complex models and solution proce- 
dures. 

In particular one can state, that in order to 
predict the critical buckling load accurately and 
to make a reliable estimate of its imperfection 
sensitivity, the nonlinear effects caused by the 
edge restraint conditions must be included in 
the analysis. Any solution procedure which 
fails to account for these effects, should be 
suspect of having provided incorrect results. 

The most approximate of the here 
described analyses, the Level-1 solutions 
which neglect the effects caused by the edge 
restraints, can still be used to great advantage 
to establish the approximate behavior of a shell 
subjected to the specified external loading. 
However, depending on the value of the 
prebuckling stiffness, resulting from the 
different types of wall constructions used, the 
solutions may be either conservative or 
nonconservative. 

As can be seen from the results shown in 
Table 5, the buckling load of the isotropic shell 
A-8 is sensitive to all the initial imperfection 
shapes investigated. For a more specific 
prediction of the final collapse load, the final 
goal of a “High Fidelity Buckling Load 
Analysis”, one has to carry out a refined Level- 
3 analysis including measured values of all the 
significant generalized imperfections such as 
the traditional shell-wall imperfections, 
variations in the shell-end or loading surface 
geometry and especially for composite shells 
variations in the shell-wall thickness 
distribution. It has been shown in Ref. 33 that 
such an approach yields very good agreement 
between the predicted collapse load and the 
experimental buckling load. 
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b Buckling Component Modes 


a. Membrane prebuckling, n=9 
=1.00006 



b. Buckling Component Modes I ! Buckling Component Modes 


b. Nonlinear prebuckling, n=24 c. Nonlinear prebuckling, n=24 

=0.844480 ^c' = 0.844481 


Buckling modes of the axially compressed isotropic shell A-8 
SS-3 B.C. -N x = -N 0 , v = W=M x =0:N<^ = -49.51517 Ib/in. 
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d. Circumferential trace of the critical (lowest) buckling mode 

Fig. 5 Buckling mode of the axially compressed isotropic shell A-8 

SS-3 B.C. N x = -N 0 , v = W = M x = 0; Nc^ = -49.51517 Ib/in; STAGS-A results. 
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d. Circumferential trace of the critical (lowest) buckling mode 


Fig. 7 Buckling mode of the axially compressed isotropic shell A-8 

SS-3 B.C. N x = -N 0 , v = W = M x = 0; Nc^ = -49.51517 Ib/in; STAGS-A results. 
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imp*rf«ct on S« r>*itlvity of liotropic SK>«II A — S — SS — 3 B.C. 

Axisymme'ric lmp«rfection: Wbar/H — xlbar t *[co*(2*pl*K/L)] 


- AXBlF.n- 9 •O-O-'O- MANJIU.ISA.rt— 9 



Fig. 8 Imperfection sensitivity for long wave axisymmetric imperfection under axial compression 
(SS-3 B.C.: N x = -N 0 , v = W = M x = 0; = -49.51517lb/in). 


Imp«rf«c ion S«sn«ftfvlty of Isotropic Sri*lt A— B — SS — 3 B.C. 
Axiiymmai ric Imp.rfsction; Wbor/h — xibar 1 -[ cos(34*pi«»x/L) ] 
AXBlF.n -27 - 0-0 * 0 - MANIL.ISA.n-iA 



Fig. 9 Imperfection sensitivity for short wave axisymmetric imperfection under axial compression 
(SS-3 B.C.: N x = -N 0 , v = W = M x = 0; N,^ = -49.51517 Ib/in). 
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Fig. 10 Contour map of the critical buckling mode for axisymmetric imperfection. 
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Aaymmatrle im p a rf ac tio n of iaotroprc A— 6 BFACT 

Wbor/h — *Jbar2-[*lnCpr-*/'l->'*coa<:®-y/«)3 



Fig. 1 1 Imperfection sensitivity for long wave asymmetric imperfection under axial compression 
using Level-1 routine BFACT (N^ = -49,51517 Ib/in). 


.rn pa '-'faction Mnaitivity of iaotropie ■ A— B uafng BFACT 

woor/h — Kit>or2»C«in( 1 7»pi«*f/L)*coaC^ 7 *y/ R 3 ^ 



Fig. 12 Imperfection sensitivity for short wave asymmetric imperfection under axial compression 
using Level-1 routine BFACT (N c ^ = -49.51517 Ib/in). 
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A»ymm«trle lmparf«etion Ivl ty of isotropic ahall A — S — 3S — ■ 3 ■ .C . 


Wbor/h — sibo •[•In (p l®w/L)*co*C8»y/W)3 

ANILISA,n-9 -0*0—0* COlLAPSC,n-C 



Fig. 13 Imperfection sensitivity for long wave asymmetric imperfection under axial compression 
(SS-3 B.C.: N x = -N 0 , v = W = M x = 0; = -49.51517lb/in). 


Aaymmatrlc tmparfactlon S*nslt iv-i ty of Isotropic A— C — SS — J S.C. 


Wbqr/h — ¥lbor2»[«ln<1 7*pl*»/L)*e»«<2*«y/R)] 

- ANIUSA.n-2* -0*0— O* Cai_l_A*SC.n— X* 



Fig. 14 Imperfection sensitivity for short wave asymmetric imperfection under axial compression 
(SS-3 B.C.: N x = -Nq, v = W = M X = 0; = -49.5151 7 Ib/in). 
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Fig. 1 5 Measured initial imperfections of Caltech isotropic shell A-8 [1]. 
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Fig. 1 6 Caltech test set-up and corresponding computer model. 
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